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Abstract 

We study in some detail the spectral phenomenology of models in which supersymmetry 
is dynamically broken and transmitted to the super symmetric partners of the quarks, leptons 
and gauge bosons, and the Higgs bosons themselves, via the usual gauge interactions. We 
elucidate the parameter space of what we consider to be the minimal model, and explore 
the regions which give rise to consistent radiative electroweak symmetry breaking. We in- 
clude the weak-scale threshold corrections, and show how they considerably reduce the scale 
dependence of the results. We examine the sensitivity of our results to unknown higher- 
order messenger-sector corrections. We compute the superpartner spectrum across the entire 
parameter space, and compare it to that of the minimal supergravity-inspired model. We 
delineate the regions where the lightest neutralino or tau slepton is the next-to-lightest super- 
symmetric particle, and compute the lifetime and branching ratios of the NLSP. In contrast 
to the minimal supergravity-inspired model, we find that the lightest neutralino can have a 
large Higgsino component, of order 50%. Nevertheless, the neutralino branching fraction to 
the gravitino and the light Higgs boson remains small, <^10~ 4 , so the observation of such a 
decay would point to a non-minimal Higgs sector. 
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1 Introduction 



Most studies of supersymmetric phenomenology have focused on models in which su- 
persymmetry is broken in a hidden sector at a scale of order Mgusy — 10 10 GeV. 
In these models, supersymmetry breaking is communicated to the visible sector by 
gravitational interactions. Such models give rise to superpartner masses in the 100 
- 1000 GeV range, and to a gravitino with approximately the same mass. The cou- 
plings of the gravitino to matter are negligible, suppressed by (E/Ms\jsy) 2 — 10~ 16 . 
The now-standard "minimal supergravity" model is of this type, with the familiar 
unification-scale boundary conditions on M , M1/2, and A . 

Recently, there has been a resurgence of interest in models where supersymmetry 
is broken at a scale of order M SU sy^10 5 GeV |] |12| . In these models, supersym- 
metry breaking is transmitted to the superpartners of the quarks, leptons, and gauge 
bosons (and to the Higgs bosons themselves) via the usual SU(3) x SU(2) x U(l) 
gauge interactions. Because gauge interactions are flavor diagonal, these models natu- 
rally suppress the flavor- changing neutral currents associated with the soft squark and 
slepton masses. 

Models with gauge-mediated symmetry breaking also have superpartner masses 
in the 100 - 1000 GeV range. However, because of the low scale of supersymmetry 
breaking, the gravitino is essentially massless; its couplings to the superpartners are 
suppressed by (i?/MsusY) 2 ^10 _6 . Typically, superparticles decay by cascading down 
to the next-to-lightest supersymmetric particle (NLSP), which in turn decays to its 
partner and the gravitino. Such decays give rise to characteristic experimental signa- 
tures, e.g. final states containing two photons and missing energy p], ||. 

In this paper we take a detailed look at the low-energy spectrum of the simplest 
models with gauge-mediated supersymmetry breaking. We begin by describing the 
models, elaborating on the issues associated with electroweak symmetry breaking. We 
use two-loop renormalization group equations and the full one-loop threshold correc- 
tions to determine the parameter space which gives rise to consistent radiative elec- 
troweak symmetry breaking. We find the Higgs boson and superpartner masses to the 
same level of accuracy. We comment on how our results are affected by higher-order 
messenger-sector corrections. We compare the spectra of the simplest gauge-mediated 
models with those from the minimal supergravity-inspired model. Finally, we identify 
the next-to-lightest supersymmetric particle and compute its lifetime and branching 
ratios across the allowed parameter space. 



2 Electroweak symmetry breaking 

In models of gauge-mediated supersymmetry breaking, the SU(3) x SU(2) x U(l) 
gauge interactions of "messenger" fields communicate supersymmetry breaking from 
a hidden sector to the fields of the visible world. In the simplest models fl[], the 
messenger sector contains a set of vector-like fields, Mj and M i; coupled to a standard- 
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model singlet, S, through the superpotential interaction 

W messenger = \SMi~Mi. (1) 

The lowest and F-components of the singlet superfield S acquire vevs through their 
interactions with the hidden fields. This breaks supersymmetry and i?-symmetry. 

To maintain the near unification of the standard-model gauge couplings, we will 
take the fields and Mj to transform in complete SU(5) representations. In the 
same spirit, we will also require that the gauge couplings remain perturbative up to 
the unification scale. This implies that we can consider at most four 5 + 5 pairs, or 
one 10 + 10 and one 5 + 5 pair. 

In what follows we will take the fields Mj and Mi to lie in (n^, nio) 5 + 5 and 10 + 10 
SU(5) representations. We will assume that these fields couple to the singlet S through 
a single Yukawa coupling A at the unification scale. (See ref. f|, for a discussion of 
variations.) Below that scale, the SU(5) representations split apart, and the Yukawa 
couplings evolve according to their own renormalization group equations 0. We will 
ignore this splitting in most of what follows; we will remark on it briefly in sect. 3. 

The lowest and F components of the superfield S acquire vevs through interactions 
with the sector which dynamically breaks supersymmetry. These vevs induce masses 
and mixings for the messenger fields. The messenger fermions gain mass M = \{S), 
while the messenger scalars obtain the mass matrix 

M 2 X(F S ) 
MF S ) M 2 

From this point on, we omit the brackets which denote the vevs. 

Following the philosophy of dynamical models, we will assume that the standard- 
model //-term and all soft supersymmetry-breaking terms arise dynamically. For the 
case at hand, the messenger fields transmit the supersymmetry breaking to the visi- 
ble sector through loop diagrams which contain insertions of the S superfield. Such 
diagrams induce weak-scale masses for the gauginos and scalars of the minimal su- 
persymmetric standard model. However, they cannot give sizable values for the soft 
supersymmetry-breaking A parameters. In what follows, we shall set all A-terms to 
zero at the messenger scale. 

The messenger fields induce the following gaugino 

M,(M) = (n 5 + 3n w )g (A) ^Al A ( 3 ) 

and scalar 

m 2 (M) = 2(n 5 + 3n 10 )/ (A) £ hCl {^-)\ 2 (4) 

masses at the scale M, where A = F s /S and ki = 1,1,3/5 for SU(3), SU(2), and 
U(l), respectively. The Cj are zero for gauge sing lets, and 4/3, 3/4, and Y 2 for the 
fundamental representations of SU(3), SU(2), and U(l). Here Y = Q — I 3 denotes the 
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usual hypercharge and we use the grand unification normalization for a\. Because a 
pair of 10 + 10 fields contributes to the soft masses as if = 3, we will set nio = 
and only consider changes in n 5 . 

The messenger-scale threshold functions || [7[ 



1 + x 

g(x) = — ^— log(l + x) + (x — >• —x) , (5) 



x 2 



X 2 



x \ ) / 2x 
log(l + x) - 2Li 2 ( — — + -Li 2 



-- , , + (x — > —x) , (6) 

.1 + xj 2 + v ; v ; 

have the property that g(x),f(x) — > 1 as x — > 0, or A <C M. In this limit, the 
expressions (^) and (|4]) take the characteristic simple forms [|| that are often associated 
with gauge-mediated models. 

The region A — > M corresponds to x — > 1, where g(l) ~ 1.4 and /(l) ~ 0.7. We 
must exclude the limit M = A because it gives rise to a massless messenger scalar. For 
the purposes of this paper, we restrict ourselves to the region x < 0.97, or M/A > 1.03. 
This corresponds to an upper limit on the fine tuning of the messenger masses; it is 
obtained by requiring that the average scalar mass-squared be less than 30 times the 
light scalar mass-squared. 

Equations (f|) and @ serve as boundary conditions for the renormalization group 
equations at the messenger scale, M. They give rise to rather generic predictions for the 
soft supersymmetry-breaking gaugino and scalar masses. In contrast, the boundary 
conditions for B(M) and yu(M) are more model-dependent .[] New interactions are 
necessary to induce these terms because they violate a Peccei-Quinn symmetry in the 
effective action [§, |9j. The new interactions can give rise to additional contributions 
to the scalar masses beyond those in (|J). In particular, they can give additional 
contributions to the Higgs masses, rn? H and rn 2 H2 . 

In this paper, we will not commit ourselves to a particular model for B(M) and 
/i(M). Instead, we will take a more phenomenological approach and treat them as 
free parameters. We will, however, assume that the soft Higgs masses m 2 Hl and m# 2 
are given by eq. (^). We will also require that electroweak symmetry be radiatively 
broken. 

Our approach is as follows. We start at the messenger scale, M, and fix the bound- 
ary conditions (|3]) and (f|). We use the two- loop renormalization group equations to 
run the soft masses down to the squark mass scale, Mq ~ y^A/QO, where we impose 
electroweak symmetry breaking and calculate the supersymmetric mass spectrum. At 
the squark scale, we consistently include all one-loop weak-scale threshold corrections. 
These corrections play an important role in determining B(Mg) and fi(Mg). We then 
run B, fi and the gauge and Yukawa couplings back to the messenger scale. We repeat 
the procedure until we find a self-consistent solution for B(M) and /i(M) in terms of 
the Z-boson pole mass, M z , and the ratio of Higgs vacuum expectation values, tan/?. 



1 The parameter \i is the standard supersymmetric Higgsino mass; B is the dimension-two soft mass that is often 
denoted -B/i in supergravity-inspired models. 
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(We take the top- and bottom-quark pole masses to be m t = 175 GeV, m b = 4.9 GeV, 
and the MS value of the strong coupling to be a s {Mz) = 0.118.) 

We note that there are messenger-scale threshold corrections which we do not take 
into account, even though they are formally of the same order as the weak-scale thresh- 
old corrections. We choose to ignore these corrections because they are model depen- 
dent. We will, however, estimate their importance by determining the sensitivity of 
our results to small changes in the messenger-scale boundary conditions. 

With these assumptions, the parameter space of this minimal model is described 
by tan/?, the scale A, the messenger scale M, along with n^, the effective number of 
5 + 5 messenger fields, and the sign of fi. We will first set = 1 and M = 2A. This 
leaves a two-dimensional parameter space, (A,tan/3), for each sign of \x. 

In Fig. |l] we show our results for B(M) / /x(M) and /x(M) in the (A, tan (3) plane. We 
see that consistent electroweak symmetry breaking requires |/z(M)| to be in the range 
from 150 GeV for small A to over 1 TeV for large A. We also find that B(M)/ii(M) 
ranges from near zero, at the largest values of tan /3, to >500 GeV at small tan f3. We 
do not consider A > 300 TeV because of fine tuning considerations. 

Various constraints exclude the border regions in Fig. [I] (and the remaining contour 
plots). The region of small tan f3 is excluded because the top Yukawa coupling diverges 
below the unification scale. This is not necessarily fatal, but since we wish to preserve 
perturbative grand unification, we require A^Mqut) < 3.5. The region of very large 
tan/3 is excluded because electroweak symmetry is not broken. Typically, we find 
m\ < in the large tan /3 region. 

The region at small A is excluded because of the usual experimental bounds on 
the masses of super symmetric particles, most notably > 65 GeV. (We use the 
bounds listed in Ref. [0.) These bounds hold when the LSP is the lightest neutralino 
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Figure 2: The parameters fi(M) and B(M) / /x(M) as a function of the scale, Q, 
where the renormalization group evolution is stopped and electroweak symme- 
try breaking is imposed. The dotted (dot-dashed) lines include the one-loop 
(two-loop) evolution equations, but no weak-scale thresholds. The dashed 
(solid) lines include the one-loop (two-loop) evolution equations, plus the full 
set of one-loop weak-scale thresholds. The dotted vertical lines on each plot 
denote the squark scales Mq/2, Mg and 2Mg. 



and is stable. They also hold in gauge-mediated models when the NLSP is Xi an d 
it decays outside the detector. In sect. 4 we delineate the regions of the parameter 
space where this applies. If the NLSP decays inside the detector, the collider bounds 
are stronger because the decay modes (e.g. X1X1 ~~ * 77+ missing energy) have much 
smaller backgrounds f2|. 

In Fig. H we plot the scale dependence of B(M)/ [i(M) and fJ,(M) to illustrate the 
importance of the weak-scale threshold corrections. We see that our results, which 
consistently incorporate the weak-scale thresholds, have significantly less scale depen- 
dence than they would if the corrections were ignored. Furthermore, we see that the 
threshold corrections are especially important at large tan (3, where B(M) is small. 
The largest threshold corrections arise from squark loops, so the threshold corrections 
are smallest in the vicinity of the squark scale, Mg = y/nsA/90. Indeed, in Fig. | we 
see that the tree- level result is nearly equal to the full result at the scale Mq/2. If we 
vary the scale over a reasonable range, e.g. Mq/2 to 2Mq, we find that the weak-scale 
threshold corrections (between 5% and 15%) are generally larger than the corrections 
from the two-loop evolution (about 2%), as might be expected because of the relatively 
small amount of running. 

The region B(M) ~ is of considerable phenomenological interest [|S], [LTJ. For 
instance, models where B(M) ~ give rise to large tan/5 without fine tuning M. 
Furthermore, since the A-terms are also small at the scale M, such models have a 
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naturally small neutron electric dipole moment [[J. Indeed, only the small 0{aM2/n) 
starting values of A(M) and B(M)/ fi(M) can contribute to the neutron EDM. The 
pieces of A and B which are generated by the renormalization group, Arc and Brg-, 
do not contribute because they are proportional to M 2 and nM 2 , so the invariant 
CP-violating phases |C1| axg(A* RG M 2 ) and ajg(B* RG ^M 2 ) vanish. 

The region of B(M) ~ is associated with the region of large tan /3. At tree level, 
we have 

B = m k ~ m k _ M 2 Z 

tan /3 — cot (3 tan (3 + cot (3 ' 

which implies that B decreases as tan /3 increases. From Fig. |l| we find, for n 5 = 1 and 
M = 2A, that B(M) is small (but non-zero) along the excluded region at large tan/3. 
(Our results should be contrasted with those of ref. [|[0|1 , which do not include one-loop 



thresholds. The authors of ref. [ 10 1 find B(M) = for tan/3 ~ 20 and /i < 0.) 

Until now we have restricted our attention to the simple case where M/A = 2. In 
principle, M/A can be much larger. A large hierarchy M ^> A can arise from a small 



Yukawa coupling or from loop factors [11]. An upper bound on this splitting can be 



obtained from the cosmological constraint that the gravitino relic density not overdose 
the universe. In the usual cosmological scenario, this translates into an upper bound 
on the gravitino mass of about 1 keV [HJ. (Note, however, that larger masses can 



be accommodated if the gravitino relic density is suppressed by a brief inflationary 
epoch.) If we assume that the singlet F-term is of the same order of magnitude as the 
largest F-term in the theory, we have 



VTin — — — = A 



G M P V A / Mp ' 



M\ A 2 



Setting A = 30 TeV and A = 1, we find that m G ~ 1 keV if M/A ~ 10 4 . Hence we 
will not consider values of M/A larger than 10 4 . Note that for A < 1 the upper bound 
on M/A is correspondingly reduced. 

As we increase M for a fixed A, there is more renormalization group evolution, so 
there are larger splittings between the soft Higgs masses. These, in turn, give somewhat 
larger values of /i(M). We show this in Figs. ||(a-b), where we plot contours of (J>(M) 
in the (A, tan f3) plane, for n 5 = 1 and M/A = 10 4 . We see that the values of /i(M) are 
slightly larger than those in Fig. [I]. The situation for B(M) is more subtle because the 
region B(M) ~ is sensitive to additive radiative corrections. In Figs. |3](c-d) we show 
our results for M/A = 10 4 . We find that B(M) = occurs for fi < with tan/3 ~ 20 
to 40, depending on A. For smaller M/A, the B(M) = contour in Fig. |3|(d) moves 
to the right, to larger values of tan/3. 

In Fig. [| we plot the pole mass of the lightest Higgs boson, m h , and of the CP-odd 
Higgs boson, m^, in the (A,tan/3) plane, for M/A = 2 and 10 4 . We see that m h <130 
GeV (see, however, Ref. |l2|), and that rriA is nearly always larger than 200 GeV. For 
such large values of rriA, all three heavy Higgs bosons are nearly degenerate in mass. 

The above results depend sensitively on n 5 , the number of 5 + 5 pairs. In Fig. |5| we 
show, for M = 2A, contour plots with = 2, 3 and 4. The change in the parameter 
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/x > M/A = 10 4 /J, < 




1 2 5 10 20 50 1 2 5 10 20 50 

tan/3 tan/3 



Figure 3: The same as Fig. [3], with M/A = 10 4 and n 5 = 1. The magnitude 
of fi(M) is increased relative to that of Fig. 1. 

/i(M) is dominated by the change in the scalar masses. Indeed, we find /x(M)/A oc y/n^, 
as expected from eq. @. The ratio B(M) / fi(M) obeys the same scaling at small tan /3. 
At large tan/5, however, the change in B(M) is more complicated. This is the region 
where B(M) is small, so its value depends sensitively on radiative corrections. For 
/i > 0, there is a line in parameter space where B(M) = at large A. For fi < 0, 
however, B(M) is never zero. We find B(M)/ fj,(M) > 10 GeV for n 5 = 4. 

The results presented here have varying degrees of sensitivity to the input param- 
eters. We illustrate this sensitivity by computing the changes in B(M)//jl(M) and 
fi(M). We first vary the quark masses, m& = 4.9 ± 0.5 GeV and m t = 175 ± 5 GeV, 
and find a shift A\B/u,\<10 GeV and \Au./u,\<7% in the region B(M) ~ 0. To es- 
timate the sensitivity to messenger-scale threshold corrections, we randomly vary the 
soft masses at the messenger scale by 5%. In the region B(M) ~ 0, we find shifts in 
B(M)/fj,(M) of up to 10 GeV, and changes in /i of order 5%. 

3 The spectrum 

In this section we find the masses of the supersymmetric particles that arise in the 
gauge-mediated scheme. We examine how the spectrum depends on the parameter 
space, and comment on the sensitivity to low- and high-energy threshold corrections. 
As above, we follow a self-consistent procedure. We start with the masses (|f) and (|J) 
at the messenger scale, M. We then use the two-loop renormalization group equations 
to run these masses to the squark scale, Mg = y/n^A/90. At that scale we apply the 
one-loop threshold corrections and impose electroweak symmetry breaking. We then 
calculate the superpartner masses, and run the gauge couplings back to the messenger 
scale. We iterate this procedure to determine the consistent one-loop superpartner 
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/x > fx < 




1 2 5 10 20 50 1 2 5 10 20 50 

tan/3 tan/3 



Figure 4: Figures (a) and (b) show contours of the light CP-even Higgs-boson 
mass, rrih, in GeV, while Figs, (c) and (d) show the CP-odd Higgs-boson 
mass, rriA, in GeV. We set n 5 = 1, and the solid (dashed) lines correspond to 
M/A = 2(10 4 ). 



pole masses. The weak-scale threshold corrections for all the superpartner masses are 
contained in ref. [13| . 



Figure |6] shows the spectrum in the canonical case where the messenger scale M is 
equal to 2 A. We plot the various masses against the mass of the lightest neutralino, 
Xi- The first point to note is that, for m^o>100 GeV, most of the curves are rather 
flat. This simply reflects the fact that most of the masses, including m^o, scale with A. 
The only exceptions are the light Higgs, h, whose mass is determined by radiative cor- 
rections (the correction to its mass-squared grows like In A), and the light tau slepton, 
whose mass is significantly affected by left-right mixing. 

Because of the large mixing in the tau slepton mass matrix, the light tau slepton 
mass is especially sensitive to radiative corrections. The one-loop self-energy correc- 
tion is typically less than (9(3%) (see Ref. Jl3|]). However, for large tan/3, the large 
threshold corrections to /i (see Fig. H) induce a significant shift in the light tau slep- 
ton mass. In Fig. |7|(a) we plot the light tau slepton mass versus tan f3. We show 
the one-loop and tree-level results, for various choices of the renormalization scale, Q. 
At large tan/3, the corrections can be of order 30%. In Fig. 0(b) we show the scale 
dependence of the one-loop and tree-level tau slepton masses, as well as the effect of 
the two-loop renormalization group evolution. The weak-scale threshold corrections 
significantly reduce the scale dependence of the pole mass. As expected, the two-loop 
renormalization group effects are small because of the small amount of running. 

The tau slepton mass matrix should be contrasted with those of the top and bottom 
squarks. These matrices can also have substantial mixing, but the mixing does not 
significantly change the mass eigenvalues. This implies that the masses of the light top 
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Figure 5: The same as Fig. |l] for varying n§. The dashed, solid and dot- 
dashed lines correspond to n 5 = 2,3,4, respectively. The values of fi(M) and 
B(M) / fi(M) scale like \/n 5 for a fixed value of A, except for B(M) ~ 0, where 
radiative corrections are important. 



and bottom squarks cannot be much less than the masses of the other squarks. Indeed, 
we find that the mass of the light top squark is typically 10 to 20% less than that of 
the other squarks, although it can be as much as 35% less at the smallest values of 
tan/3. Similarly, the light bottom squark mass is usually to 10% less than that of 
the other squarks, but it can be as much as 25% less at the largest values of tan (3. 

In Fig. |] we see that for m^o<100 GeV, the curves leave their linear trajectories, 
especially for tan/3 = 2 with \i negative. This is because for small m^o, the mixing in 
the neutralino mass matrix become increasingly important, so the Xi mass no longer 
scales with A. This is illustrated in Fig. |8], where we show the Xi mass as a function 
of A. 

The previous spectra were all computed with M = 2 A. If M is increased with 
respect to A, the masses change in two ways. First, the initial conditions are different 
because the masses at the scale M depend on ojj(M). Second, more running is needed 
to reach the squark scale, Mq. 

These effects almost cancel for the gaugino masses because the gaugino masses 
obey the same one-loop renormalization group equations as the gauge couplings. They 
do, however, change the scalar masses. For squarks, the two effects work in opposite 
directions. It turns out that for 715 = 1 the change in the boundary conditions is 
more important, so the squark masses are reduced. For = 3 the two effects largely 
cancel. For sleptons, the effects work in the same direction, and give a slight increase 
in the slepton masses. The shifts in the scalar masses are illustrated in Fig. || The 
small changes indicate that it is safe to ignore any factor of two or three change in the 
messenger quarks scale with respect to the messenger lepton scale, such as might be 



9 




40 100 200 400 40 100 200 400 40 100 200 400 40 100 200 400 
m~c (GeV) m~o (GeV) m~o (GeV) m~o (GeV) 



Figure 6: The ratio of various superpartner masses to rn^o, versus rn^o, for 
M = 2 A and n.5 = 1. The different lines for each particle correspond to the 
cases of small/large tan/3 and positive/negative \i. Except for uia and m T1 , 
the large tan (3 curves are essentially independent of the sign of /i. 

induced by renormalization group evolution of the messenger Yukawa coupling A. 

Figure [| also illustrates the importance of the threshold functions f(x) and g(x). 
In the region 1.03<M/A<2 the threshold functions lead to significant increases in the 
gaugino masses, and to small decreases in the scalar masses, relative to the case when 
M/A>2. The threshold corrections are so important that the light tau slepton becomes 
the NLSP for M near A. 

In Fig. [| we also show the effect of increasing n§ from 1 to 3. The gaugino masses 
scale like n.5, while the scalar masses go like \/n 5 . Because of this fact, for the case 
of a 10 + 10 pair of messenger fields (n 5 = 3), the light tau slepton is the NLSP over 
most of the parameter space. We discuss this in more detail in the next section. 

Note that the squark masses are light enough that they will be produced in the 
next generation of colliders over most of the parameter space. 

As above, we can estimate the sensitivity of the supersymmetric spectrum to 
messenger-scale thresholds by varying the soft masses at the messenger scale by 5%. 
We find that most of the pole masses vary by about 5%, except for the light Higgs, 
which varies by <1%, and the top and bottom squarks, which change by up to 10%. 
The light tau slepton mass varies by up to 15%, except at a few exceptional points, 
where the variation can be as large as 70%. This large sensitivity in the tau slepton 
mass occurs at large tan /3, in the region where it is potentially the NLSP. 

We will conclude this section by contrasting the predictions of the gauge-mediated 
models with the predictions of the minimal supergravity model. At the unification 
scale, the inputs to the supergravity model are a universal scalar mass, M , a common 
gaugino mass, M1/2, and a trilinear scalar coupling A . As in the gauge-mediated 
model, radiative breaking of electroweak symmetry breaking is imposed. 
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Figure 7: (a) The light tau slepton mass versus tan/?, for 715 = 1, A = 100 
TeV, M = 2A, and \i < 0. The one-loop result (solid) and tree-level results 
for Q = Mg/2, Mq, and 2Mq are shown. The tree-level curve for Q = Mg/2 
is closest to the full result, (b) The scale dependence of the tau slepton mass 
for the same parameters, with tan (3 = 50. The full result (solid) is contrasted 
with the tree-level result (dot-dashed). The corresponding curves obtained 
using one- loop renormalization group evolution are shown (dashed, dotted). 



In general, the two models predict qualitatively different spectra. The supergravity 
model includes separate scales for the gauginos and the scalars, so the scalars can be 
much heavier than the gauginos. The minimal gauge-mediated model has only one 
scale, so the ratios of the gaugino to the scalar masses are more or less fixed. We 
illustrate this in Fig. |T0|(a), where we show the ratio of the (first or second generation) 
squark mass to the gluino mass in the two models.0 In the supergravity case, this 
ratio varies from about 0.9 to 5. For gauge-mediated models, the ratio depends on n 5 . 
For n.5 = 1, the squarks are 20 to 40% heavier than the gluinos, while for = 4 the 
squarks are about 5 to 10% lighter. (These numbers and the bands on the figure ignore 
the region M<2A, where the ratio mq/m g falls by 20%.) Note that gauge-mediated 
models with = 4 predict the same ratio of squark to gluino mass as the supergravity 
models with M <Mi/ 2 . 

In the region where the two models predict the same gluino and squark masses, they 
do not typically predict the same slepton masses. We illustrate this in Fig. 0(b), where 
we show a supergravity scatter plot of the ratio of the left-handed to the right-handed 
selectron mass, versus M /M 1 / 2 - On the scatter plot we superimpose the prediction 
of the gauge-mediated model, for n 5 = 4, divided up into regions of small and large 
right-handed selectron mass (m^ < M z and m e - fl > M z ). In the small region, the 
D-term contributions to the slepton masses are important. The D-term contributions 
can make the selectron ratios coincide with those from the supergravity model in the 
region M Q <C M]/ 2 if tan/3 is less than 3 in supergravity model and larger than 5 in 
the gauge-mediated model. 

2 The supergravity scatter plots were generated as in Ref. B . 
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The ratio m^o/A versus A, for the same choice of parameters as in 



If we go even further, we find that the two models can predict the same gluino, 
squark, charged slepton and light chargino masses. We give an example in Table 1. 
This table lists the spectra for the supergravity model parameters tan (3 = 2.1, Mq = 5 
GeV, Mi/2 = 145 GeV, A = 686 GeV, and /x < 0, and for the gauge- mediated model 
parameters tan/5 = 14, A = 11 TeV, M/A = 3716, 7x5 = 4 and /x < 0. The overlap 
can occur only for small selectron and Xi masses, where the slepton .D-terms and the 
gaugino/Higgsino mixing come into play. 

Table 1 also illustrates that the models differ in their predictions for other observ- 
ables. If we scan over the parameter spaces of the two models and match the gluino 
masses to 8%, and the light chargino and (first and second generation) squark and 
selectron masses to 2%, we find that the two models predict different values of tan/3, 
as well as significantly different masses for the neutralinos, the sneutrinos, the heavy 
chargino and the light Higgs boson. They also predict different values for the Higgs 
parameters /x and tua, and, typically, different masses for the third-generation squarks 
and sleptons. 

We illustrate these differences in Fig. |TT|, where we show scatter plots of the values 
of tan (3, rrih,m a, and m~+ in the two models. We see that the models are clearly 
distinguished. Note that since both models have essentially one scale, there is a strong 
correlation between the heavy Higgs and Higgsino masses. (The two scales are A 
and Mi/2, since M — 0.) We conclude that it is possible, but not probable, that 
measurements of the gluino, squark, slepton and light chargino masses may not be 
enough to discriminate the gauge-mediated from the supergravity model. However, 
additional measurements would easily rule out one (or both) of the models. 
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Figure 9: Various sparticle masses versus the ratio M/A, for A = 50 TeV, 
tan/3 = 2 and \i > 0. The solid (dashed) lines correspond to n 5 = 1 (715 = 3). 



4 Properties of the NLSP 

In the preceding section, we saw that the NLSP is either the light tau slepton, fx, or 
the lightest neutralino, Xx- Over most of the parameter space, the number of 5 + 5 
pairs determines the NLSP. For n.5 = 1, the Xx tends to be the NLSP. For larger 715, 
the balance tips towards the fi because as n 5 increases, the scalar masses increase less 
than their gaugino counterparts. 



In Fig. [12] we plot contours of m^o = in the (m^o, tan/3) plane, for each value 
of 715 between 1 and 4. Each subfigure contains four contours, two for each sign of /i, 
and two for M/A = 2 and 10 4 . We see that for = 1, the neutralino is the NLSP, 
except for small regions at large tan/3>25 and small Tn^o. For n 5 = 2 and M/A = 10 4 , 
X? is the NLSP for tan/?<20 - 30. For n 5 = 2 and M/A = 2, the light tau slepton 
is the NLSP over most of the parameter space, as it is for n§ = 3 or 4 and arbitrary 
M/A. For ri5 > 2, the neutralino is the NLSP at very small m^o<100 GeV. 

These results are subject to significant uncertainties from the unknown messenger- 
scale thresholds. In particular, the line where m^o = is sensitive to potentially large 
one-loop corrections to the mass of the light tau slepton. We find that 5% variations 
in the messenger-scale boundary conditions give rise to a m^o — m fl mass difference of 
up to 30 GeV in the region where the difference is less than 100 GeV. 

The collider phenomenology in gauge-mediated models depends crucially on the 
nature of the NLSP, especially when it decays inside the detector. The lifetime of the 
NLSP depends on its mass, m, and on the ratio M/A, as follows, 



FS 1 (M 
t oc — -a 1 . . , 
m° m VAA/ 



:) • < 9 > 

For the case at hand, where the NLSP is either the B-ino or the fx, we can put in the 
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Figure 10: (a) The ratio of the (first or second generation) squark mass to the 
gluino mass in the supergravity-inspired (dots) and gauge-mediated models 
(bands). The four bands in the gauge- mediated case correspond to n 5 = 
1,2,3,4. The largest ratio (uppermost band) is for n 5 = 1. (b) The ratio of 
the left-handed to the right-handed selectron mass in the two models. The 
scatter plot shows the supergravity model prediction. The bands mark, in the 
gauge-mediated model with n 5 = 4, the regions where mg R < Mz (grey) and 
m e - fi > M z (hatched). 



appropriate factors of it and ai and write the lifetime in units of meters, 

/100 GeV\ f M\ 2 - 
cr ~ f — J f — J x 10" 5 meters. (10) 

For = yj ' E 2 /m? — 1 ~ 1, this is the approximate decay length. We see that for 
fij ~ 1, A ~ 1 and M/A<100, the NLSP will decay inside the detector. However, for 
sufficiently small A it will decay outside the detector. 

In Fig. [H| we illustrate the neutralino lifetime in different regions of parameter 
space. In Fig. |l3"|(a) we plot the lifetime versus the neutralino mass, for M/A = 2 and 
100, with 77,5 = 1, A = 1 and tan/5 = 2. In Fig. 0(b) we show a scatter plot of the 
neutralino lifetime versus M/A. The lifetime increases by eight orders of magnitude 
as we increase M/A from 1 to 10 4 . In this and the next scatter plot, we vary the 
parameters n§ from 1 to 4, A from 16 TeV/n.5 to 300 TeV/ v/ns, M/A from 1.03 to 
10 4 , and tan/? from 1.2 to 70. The latter three variables are sampled on logarithmic 
measures, subject to the phenomenological constraints discussed in sect. 2. 

If the tau slepton is the NLSP, it will decay to a r lepton and a gravitino with a 
branching fraction that is essentially 100%. The lightest neutralino has many more 
decay modes because it has Higgsino and gaugino components. The x\ can decay to 
either a Higgs or a gauge boson, plus a gravitino, G. In Fig. [14] we show the Higgsino 
and photino components of the Xi- 111 Figs. |l4|(a) and |l4"l(b) we set n 5 = 1 and see 
that the NLSP is approximately 60 to 90% (40% to 70%) photino for fj, > (p < 0). 
The Higgsino component is quite small. The Higgsino component can be larger if n 5 
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Table 1: Comparison of the spectra of a gauge- mediated (GMSB) and 
supergravity-inspired (SUGRA) model. The masses are in units of GeV. 



is larger, because is reduced relative to the gaugino masses. In Fig. |14](c) we see 
that for 715 = 4 and /i > 0, the Higgsino component of the NLSP can be large. In fact, 
it can be as large as 50%. 

The x° partial widths can be readily computed. We express them in units of 



A 



96ttM?M| 



m% 
32ttF| ' 



(11) 



We first write the well-known decay rate to a photon and the gravitino, 

A^Tffi - G 7 ) = 2k it (12) 

For the remaining decay modes, we give the three body decay formula to account for 
the decay to standard model fermions via real or virtual boson exchange. If we sum 
over final-state standard-model fermions, we find 



A -i T ($ -> Gf+f- 



(Z) + K iZ L Iq (Z) 



(13) 
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Figure 11: Scatter plots of observables in the regions of parameter space where 
the gauge-mediated and supergravity models predict the same gluino, squark, 
slepton and light chargino masses (see text). Figure (a) shows a scatter plot 
of rrih versus tan (5 in the gauge-mediated (GMSB) and supergravity models 
(SUGRA). Figure (b) shows m~+ versus m^. 



where the sum I]/ is over all fermions with mass less than m^o/2, and N£ is the number 

of colors: 3 for quarks and 1 for leptons. In the virtual photon contribution, s min is a 
detector dependent cut-off on the fermion pair invariant mass. Also, we define 

(14) 





\N a c w + N i2 s w \ 2 , 




KiZ T = 


\N i2 c w - Nns w 2 , 




KiZj = 


c w s w (\N i2 \ 2 - \N U 


I 2 ) + (c 2 w - *w)Ke (TV* N i2 ) , 


K iZ L = 


\N i3 C/3 - Ni 4 S/3\ 2 , 




^iA = 


\N i4 Cf3 + N i3 S/3\ 2 , 




KiH = 


\N i3 c a + N i4 s a \ 2 , 




K ih = 


\NnC a — N i3 s a 2 , 




KiHh = 


SaC a (\N i4 \ 2 - \N i3 f 


) + (c 2 a - sDUe (N* 3 N i4 ) . 



In writing the H-h interference contribution, we have neglected all Yukawa couplings 
except those of the the bottom and tau. The 7-Z {H-h) interference terms go to zero 
in the limit m^o ^> Mz (m^o ^> m#). Note that if m^o > M v (if = h, H, A), the decay 

rate A~ l T(x°i — ► Gip) = K iip Ii(ip). 

In these expressions, the integrals I n are given by 

with = M 2 /m? and e v = r^M^/mL These integrals reduce to (1 — R^) 4 for 



16 



20 25 30 35 40 45 1 2 5 10 20 50 



200 




(o) n 5 = 3 






(d) n 5 = 4 


100 




%v / 
/ v / 

/'■'' Y° ~ ~ - V ' / 






\ " ' X 1 ~~ v< . / 


50 




\ 1 / 














[ i i i 1 1 1 1 i , 


1 


2 5 10 20 50 

tan/3 


l 


2 5 10 20 50 

tan/? 



Figure 12: Contours marking the boundary where the mass of the light tau 
slepton is equal to the mass of the lightest neutralino, for n.5 = 1,2,3 and 4. 
The solid (dashed) line corresponds to /J > 0, M/A = 2 (10 4 ). The dot-dashed 
(dotted) line corresponds to /x < 0, M/A = 2 (10 4 ). The excluded region is 
plotted for /j > and M/A = 2. The label x? (fi) marks the region where the 
neutralino (tau slepton) is the NLSP. 



small (i.e. in the narrow width approximation). The integral Iio(h,H) is 

[x - R H ) (x - R h ) + e h e H 



ho(h,H) 




x) 4 



((x-R h y + el)((x-R H y + el) 



(16) 



Our results complete the formulae for the two-body decay rates of the Xi given in 
ref. ||. Our results include the contribution of the virtual photon and the contributions 
from Z-7 and H-h interference. Note that the formula for T(x° — > G+ Higgs boson) 
contains the function Ii, while the formula for T(x° — > G+ longitudinal Z boson) 
contains the function Iq. At first sight one might think that this violates the electroweak 
equivalence theorem, since the decay rate to the longitudinal Z is not equal to the decay 
rate to the Goldstone boson. However, the equivalence principle does indeed hold in 
the applicable regime, since the two functions approach each other for m^o 3> Mz- 

Xi 

In Fig. |15|(a) we show the branching fractions of the Xi versus m^o, for tan/5 = 2, 
= 1, and M = 2A. (We plot the magnitude of the Z-7 interference contribution, 
|r^ 7 |/r tot , and we take into account 5 quark flavors.) The decay to the photon domi- 
nates, as expected, since the Xi has a large photino component. The branching fraction 
to fermions via an off-shell photon varies from about 3 to 4% for s^ in = 1 GeV. 

In the parameter space associated with the minimal models, the branching ratio to 
fermion pairs via heavy Higgs boson exchange is negligible. For the parameters corre- 



sponding to Fig. 15, the heavy Higgs exchange branching ratio is of order 10 . The 



branching ratio associated with H-h interference and the branching fraction associated 
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Figure 13: (a) The lifetime of the NLSP, for tan/3 = 2 and n.5 = 1, with 
/i > (solid line) or /1 < (dot-dashed). For fixed r« x o, the lifetime scales 
like (M/A) 2 . (b) A scatter plot of the neutralino lifetime, restricted to cases 
where the neutralino is the NLSP. 



with virtual /i-boson exchange are about an order of magnitude larger. 

Because of the kinematical suppression, the branching fraction to the Z-boson rises 
slowly above threshold, to about 15% (20%) for m^o = 2M Z (3M Z ). We illustrate this 
in Fig. |15|(b), where we show a scatter plot of the branching fraction of the lightest 
neutralino to a gravitino and a fermion pair. We see that the branching ratio depends 
only on the neutralino mass. This would not be true if the part of parameter space in 
which the Xi has a large Higgsino component (715 = 4, // > 0, small A, see Fig. 0(c)) led 
to a sizable branching ratio to the Higgs boson. However, in this region of parameter 
space either the f\ is the NLSP or the xi mass is below the /i-boson threshold. Hence, 
the associated partial width is negligible. 

As discussed above, the mechanism that generates B(M) and n{M) might result 
in extra terms to the scalar masses beyond those in eq. (§). In this case, the NLSP 
can decay predominately into a Higgs boson. To study this possibility, we take /1 and 
B to be independent low-energy parameters, and focus our attention on the case of 
small fi, so the lightest neutralino is predominantly Higgsino. In Fig. 16 we plot the 
lifetime and branching fractions of the Xi as a function of m^o , for various choices of 
parameters. Note that in the region m^o<250 GeV the lightest neutralino is largely 
Higgsino, and m^o — 

The lifetime and branching fractions of the NLSP in the Higgsino region are shown 
in Fig. 16. In Fig. [16|(a) we see that the lifetime varies by 6 or 7 orders of magnitude as 
m^o varies over less than one order of magnitude. This can be contrasted with the case 
where the neutralino is dominantly I?-ino, where the lifetime varies by less than one 
order of magnitude (see eq. @ and Fig. 13). The branching fractions of the lightest 
neutralino are shown in Figs. |16|(b-d). 

These results can be understood as follows. In the small fi region the NLSP is 
primarily Higgsino, so it prefers to decay to a Higgs boson. Below h and Z threshold, 
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Figure 14: Figures (a) and (b) [(c) and (d)] show the photino component of 
the lightest neutralino, K\ 1 [the Higgsino component, Kih], in per cent, for 
M/A = 2 and n 5 = 1 [n 5 = 4]. 



however, it is forced to decay through its suppressed photino component. This accounts 
for the long lifetime and the large photino branching fraction for small m^o. As the 
mass increases, the h and Z channels open, so the lifetime decreases and the branching 
fractions are determined by n\z L and Kih. For small tan /3, the decay to h (Z) dominates 
for /i positive (negative). For large tan/5, K\z L — for each sign of fi. In this case 
the branching fraction to Z is larger because the Higgs is heavier, so Iq{Z) > Ii(h). 
Finally, as m^o increases still further, the decay to photons again dominates, because 
in this region, the NLSP is predominately photino. For m^o ~ 300 GeV, the decay 
length approaches 10~ 5 meters. 

The neutralino lifetime in the Higgsino region can be readily scaled for other values 
of ri5 and M/A. As before, for fixed m^o, cr scales like (M/A) 2 . For fixed M/A and 
fixed gaugino masses, it scales like 1/n 4 because r oc Fj//i 5 oc A 4 /// 5 . 

These results should be contrasted with those of the minimal models, where we 
found that, above h threshold, the branching ratio to the h varies from 10~ 8 to 
10~ 4 . Therefore collider events with four 6-jets and missing energy would imply a 
non-minimal Higgs sector in these models. 

5 Conclusions 

In this paper we examined the detailed low-energy spectrum of gauge-mediated su- 
persymmetry breaking models. We used two-loop renormalization group equations for 
the gauge and Yukawa couplings, and for the soft supersymmetry breaking parame- 
ters. We imposed consistent one-loop radiative electroweak symmetry breaking under 
the assumption that the mechanism that generates the /i and B parameters does not 
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Figure 15: (a) The branching ratios of the Xi i n the case M/A = 2, n 5 = 1 
and pi > 0. The solid line near 1 is the branching ratio to a photon and the 
gravitino, and the dashed line is the branching ratio to fermions via a virtual 
photon, (b) A scatter plot of the neutralino branching fraction to fermion 
pairs, in the cases where the neutralino is the NLSP. 



induce any extra contributions to the scalar masses. 

We examined the phenomenology in the case of an arbitrary number of 5 + 5 mes- 
senger fields, and in the case that the messenger scale M is greater than A. We began 
by examining electroweak symmetry breaking. We initially considered the minimal 
case, n 5 = 1 and M = 2A. In the parameter space we considered, with A < 300 
TeV, we found that |/i(M)| varies from 150 GeV to over 1 TeV. We also found that 
\B(M) I fJ>(M) | varies from near (but not equal to) zero to about 500 GeV. We found the 
phenomenologically interesting region B(M) ~ to be compatible with electroweak 
symmetry breaking in the region of large tan/3, for either sign of \x. For larger M/A, 
we found B(M) ~ requires /i < 0. 

We examined the spectrum and illustrated how it depends on n 5 and M/A. We 
found that the spectrum is not qualitatively affected by an increase in M/A. For 
n 5 > 1, we found that the light tau slepton is the NLSP over most of the parameter 
space. It decays to a tau lepton and a gravitino. 

By varying the boundary conditions of the soft parameters at the messenger scale 
by 5%, we determined the sensitivity of the supersymmetric spectrum to higher-order 
messenger-sector corrections. We found that in the region B(M) ~ 0, B(M) varies 
by 10 GeV. Furthermore, the spectrum varies by 5%, except for the light Higgs mass, 
which is essentially unchanged, and the third-generation squark and slepton masses. 
In the most extreme case, we found a 70% variation in the light tau slepton mass, 
which implies a substantial uncertainty in the identification of the NLSP. 

We compared the predictions for the spectra in the gauge-mediated and supergravity- 
inspired models. In general, the two models predict qualitatively different spectra. We 
found that it is possible for the two models to give the same gluino, light chargino, and 
(first and second generation) squark and slepton masses. Such a match is only possible 
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Figure 16: The (a) lifetime and (b-d) branching fractions of the NLSP, for the 
case where the Xi is primarily Higgsino, with = 1, A = 1 and A = 200 
TeV. The m^o values are obtained by varying \i from 50 to 350 GeV. The solid 
(dashed) lines are for tan/3 = 2, with fi positive (negative). The dot-dashed 
lines correspond to tan/5 = 40, /i > 0. (The curves for tan/5 = 40, \i < are 
almost identical.) 



for very light masses, where D-terms and gaugino/Higgsino mixing is important. In 
this region, the models can be distinguished by their predictions for other observables, 
such as tan (3, as well as the sneutrino, neutralino, heavy chargino, and Higgs boson 
masses. 

For n 5 = 1, we found that the Xi is usually the NLSP. The Xi decays to the gravitino 
and either an (on- or off-shell) gauge or Higgs boson. We derived the neutralino 
decay rate and branching fractions. The lightest neutralino tends to be gaugino-like, 
and decays to a photon and a gravitino. However, in the region m^o^lOO GeV, the 
branching fraction to the Z can be larger than 20%. Leptons are easier to track 
than photons, so this mode has the potential to permit a precise neutralino lifetime 
measurement. We found that for n 5 = 4 the Higgsino component of the Xi can be as 
large as 50%. Nevertheless, over the entire parameter space the branching fraction to 
the light Higgs boson is less than 10~ 4 . 

We also examined the branching fraction of the neutralino in the Higgsino region. In 
this case the branching fraction to the on-shell /i-boson can reach over 80% well above 
threshold. Since this Higgsino region can only occur for models with non-minimal 
Higgs sectors, the observation of 4 fe-jets plus missing energy would be an important 
step towards understanding the origin of the \x term. 

Models with gauge-mediated supersymmetry breaking offer the appealing possi- 
bility that the origin of supersymmetry breaking is experimentally accessible. Once 
supersymmetry is discovered, detailed study of the supersymmetric particles, along 
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the lines suggested here, might well prove to be the first step towards uncovering the 
mechanism of supersymmetry breaking. 



Acknowledgements 

One of us (D.M.P.) thanks S. Thomas and J. Wells for useful conversations. 

References 

[1] M. Dine and A. Nelson, Phy. Rev. D 47, 1277 (1993); M. Dine, A. Nelson and 
Y. Shirman, Phys. Rev. D 51, 1362 (1995); M. Dine, A. Nelson, Y. Nir and Y. 
Shirman, Phys. Rev. D 53, 2658 (1996). 

[2] D.R. Stump, M. Wiest and CP Yuan, Phys. Rev. D 54, 1936 (1996); S. Di- 
mopoulos, M. Dine, S. Raby and S. Thomas, Phys. Rev. Lett. 76, 3494 (1996); 
S. Ambrosanio, G.L. Kane, G.D. Kribs, S.P. Martin and S. Mrenna, ibid 76, 3498 
(1996); S. Dimopoulos, S. Thomas and J.D. Wells, Phys. Rev. D 54, 3283 (1996). 

[3] S. Ambrosanio, G.L. Kane, G.D. Kribs, S.P. Martin and S. Mrenna, Phys. Rev. 
D 54, 5395 (1996). 

[4] A.E. Faraggi, Phys. Lett. B 387, 775 (1996). 

[5] S.P. Martin, |hep-ph/9608"22ll , to be published in Phys. Rev. D 55. 
[6] C. Carone and H. Murayama, Phys. Rev. D 53, 1658 (1996). 



[7] S. Dimopoulos, G.F. Giudice and A. Pomarol, preprint CERN-TH/96-171, [hep- 
ph/9607225| . 



[8] G. Dvali, G.F. Giudice and A. Pomarol, Nucl. Phys. B 478, 31 (1996). 

[9] M. Dine, Y. Nir and Y. Shirman, preprint SCIPP-96-30, |iep-ph/ 96073971 to be 
published in Phys. Rev. D 55. 

[10] K. Babu, C. Kolda and F. Wilczek, Phys. Rev. Lett. 77, 3070 (1996). 

[11] T. Hotta, K.-I. Izawa and T. Yanagida, Phys. Rev. D 55, 415 (1997). 

[12] A. Riotto, O. Tornkvist and R.N. Mohapatra, Phys. Lett. B 388, 599 (1996). 

[13] D. Pierce, J. Bagger, K. Matchev and R.-J. Zhang, preprint SLAC-PUB-7180, 
JHU-TIPAC-96011, |hep-ph/96062TT| , to be published in Nucl. Phys. B. 



[14] M. Dugan, B. Grinstein and L. Hall, Nucl. Phys. B 255, 413 (1985); S. Dimopoulos 
and S. Thomas, Nucl. Phys. B 465, 23 (1996). 

[15] H. Pagels and J.R. Primack, Phys. Rev. Lett. 48, 223 (1982). 



22 



